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>■ ' We calculate the two-loop contributions to the electric dipole moments of the electron and 

in '. 

\^ . the neutron mediated by charged Higgs in a generic supersymmetric theories. The new 

Q ' contributions are originated from the potential CP violation in the trilinear couplings of 
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the charged Higgs bosons to the scalar-top or the scalar-bottom quarks. These couplings 
did not receive stringent constraints directly. We find observable effects for a sizeable 



p • portion of the parameter space related to the third generation scalar-quarks in the minimal 

Qh' supersymmetric standard model. 
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It is widely believed theoretically that supersymmetry will play a significant role at 
certain fundamental scale that can be probed in the near future. As a result, it is one of 
the major goal of the new and the future colliders to look for signature of supersymmetry. 

The strongest constraint on supersymmetry is from the flavor changing neutral current 
process such a.s K — K mixing and /i -^ 67, or from the flavor neutral CP violating 
quantities such as the electric dipole moments (EDM's) of the electron and the neutron|jl|]. 
These quantities put strong limits on new parameters related to the first two generations 
of quarks and leptons in any generic supersymmetric theories resulting in the famous 
SUSY flavor problem and the SUSY CP problem. A solution will require theory to 
explain why these parameters are so small. One scenario out of these problems is to 
assume that the first two generations of squarks and sleptons are heavy, while keeping the 
third generation relatively light in order to produce interesting physical consequence such 
as low energy electroweak baryogenesis0. Such scenario is consistent with the sparticle 
spectrum resulting from renormalization group evolution starting from an high energy 
unified theory. It is therefore interesting to look for effective constraint on the new SUSY 
parameters related to the third generation using the low energy data or collider searches. 

In a recent paper 0, it has been shown that by considering the two-loop contributions 
to the electric dipole moment of the electron and the neutron, it is possible to put non- 
trivial constraint on the CP violating parameters related to the third generation directly. 
In fact it was shown later [Q that such constraint actually helps rule out a scenario of the 
supersymmetric theory in which the CP violating e parameter in the kaon system was 
proposed to be originated purely from a new source beyond the standard model. 0. 

In Ref.0, only the two-loop diagrams with the second loop mediated by a neutral 
Higgs boson have been considered. These diagrams are of Barr-Zee type@]. However, 
there are scenarios in supersymmetric theories in which other scalar bosons, such as the 
charged Higgs boson, may be light enough to produce larger effect. In this paper we 
wish to extend the earlier analysis to include two-loop contributions with the second 
loop involving the exchange of a charged Higgs boson. It is possible to analyse the two 
sets of two-loop contributions separately because they are gauge invariant independently. 
The charged Higgs set is more technically sophisticated and has only been studied in 
non-supersymmetric theories in Ref.0. 
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Fig. 1. Structure of the two-loop EDM diagram via Barr-Zee mechanism in 
supersymmetric theories with charged Higgs exchange in the second loop. 

The two-loop diagrams are shown as in Fig. 1, in which both the stop and the sbottom 
are involved in the first shaded loop. In a generic gauge, there are also other diagrams 
involving the unphysical Higgs boson G^. Ifowever it is known|^ that in R^ nonlinear 
Landau gauge (NLLG), such diagrams will not contribute. For this reason, our calculation 
is based on NLLG. 

We shall first discuss the first loop that creates an H'jW* vertex and then integrate 
the second loop in the EDM calculation. 



The first loop of HjW* Vertex 

The relevant interaction Lagrangian in our study is given by, 

C = vH-*i*£b + ^2tan/?m, _^^^^_, _ ^^-*^^-^^.^^ ^ i*ICtH) + H.c. (1) 

plus mass terms. Here t and b carry indices 1, 2 labeling the corresponding mass eigen- 
states which are superpositions of the bosonic chirality states, 

t)="-(t) ■(;:)=«(:-:)■ 

These unitary transformation Ut and Ub diagonalize the mass matrices 



Ml +ml + Ag^ Ml cos 2/3 - {fngii + m^ A* 

-{m,li* + m^Aq) M^ + m^ + A,^M| cos 2(3 






with q = t,b with rhi, = mfetan/5, rht = nit cot P, ^q = T^ — eq?,iv? 6y/ and the SU{2) 
relation Mi^ = M-^ . The complex phase can be factored out as 



c,= : ' ' . (4) 




with Cq = COS 6q, Sg = siii 6q for q = b,t, and 

TJig^* + iriqAq = Irriqfi* + mqAqle'^" . (5) 

The charged current couphng matrix /C in this basis of the i-b sector is real, 

V y y stCb stSb J 

The tri-hnear scalar terms for charged Higgs boson can be written as 



^2 . / -M^sin2/5 + mgtan/5 + m?cot/3 rubH - Alrhb . 



v^ \ rritfi* — Atrfit 2mbrat/ sin 2/3 



The imaginary entries are CP violating, 



^/2 2mbmt ( StSb StCb 

-ctSbXb + stCbXt ctCbXb + SfSbXt 
-SfSbXb - ctCbXt stCbXb - CfSbXt 



where 



K= ,^!^^ . i^KO ■ (9) 

w^ sm p cos p ^ ^ 
CP violation is a result of the mismatch in the phase between the S vertex of the charged 
Higgs boson and the /C vertex of the W boson. Our convention happens to give rise to 
real /C and complex S. Note that while there are two stop and two sbottom eigenstates, 
we actually need only one of them each to give contribution to the EDM of leptons or 
quarks. Therefore, one can in principle consider the numerically simpler limit in which 
one of the squark of each flavor is much lighter than the other one. Also, even if the squark 
mass matrices are accidentally real, the CP violation can still be generated through the 
weak-basis-charged-Higgs-coupling matrix in Eq.(|^). In the first term of Eq. (|^), the CP 
violation is due to the relative phase between stop and sbottom mass matrices, while in 
the second term CP violation is due to the relative phase between the weak-basis-charged- 
Higgs-coupling matrix and squark mass matrices. Note, however, that if Xt = Xb = 0, the 
phase of /i may not be zero, however sin((56 — 6t) will be zero and there will be no CP 
violation in the charged Higgs interaction in our approximation (of ignoring the squark 
flavor mixing). 

We work out the 1-loop amplitude for the process H{k + q) ^ ■y{k,fi) + W*{q,a). 
Besides the 3-prong irreducible diagram. Fig. 2(a), 
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Fig. 2 One loop H 
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W 7 Feynman diagrams in the non-linear Landau gauge. 



there is a one-particle reducible H~-W~ bubble diagram, Fig. 2(b), with the external 
photon k attached to the outgoing leg W~ . This amplitude is not vanishing even at the 
limit of the gauge parameter ^ ^ in our calculation based on the R^ non-linear Landau 
gauge (NLLG), where the unphysical Goldstone boson coupling in GW~'y is absent. One 
notices that the W propagator has the form 



p-/3(p) 



<T/3 , P"p'^(g-1) 



g<rp ^ 



<M^ c^o g 



o"/3 n^nl^ 



pCTpP 1^2 ^p'^pl^ 



(10) 



p2 -Mir p2 _ M^ (p2)2 

The first term gives vanishing result when it is convoluted with the bubble amplitude 
of the form -Bpo-. The second term will give finite result even when ^ — i> as there are 
singular terms from the tri-gauge-boson coupling in the Non-Linear Landau gauge. 
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where we have extracted the leading constant term in the limit ^ ^ 0. We need also 
to include a sea-gull diagram. Fig. 2(c), with the external photon k sticked to the t- 
h-W vertex. The photon (fc, yu) is attached to the external electric field for the EDM 
measurement. We only keep the lowest order of k in our calculation. 
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Note that the amplitude in NLLG satisfies the Ward identity k^T^" = even when 
q' + M^. 



The second loop and EDM 

The final 2-loop EDM amplitude (Fig. 1) of the lepton becomes 
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(13) 



Here we integrate Q^ over (0, cxd), and y over (0,1). In our EDM convention, M. 
—a^'^ku'j^di, we have the contribution to EDM from the charged Higgs sector. 
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The above integral is derived based on the relation. 
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(15) 
(16) 



(x + a) (x + 6) (x + c) a — h 

The value of JFg is negative and approaches zero as the arguments r, (3 approach infinity. 
Therefore JF is also nagetive. Note that Y.c^d^'^i.^l d)^c,d = as expected because if 
Mq^ = Mq^{q = t, b) the mass matrix A^? is proportional to the diagonal matrix and the 
CP violating phase 6q disappears. In that case, even though the charged Higgs coupling 
matrix in Eq.(0) may appear to be complex, but the off diagonal phases can easily be 
absorbed into tn and 6r. Note also that if M^^ = M^^, the CP violation will depends on 
At only and the contribution from the first term in Eq.(|D will disappear as expected. 



Comparison with Pseudoscalar Contribution 

The neutral A^ couples to stop and sbottom in the form 
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The real part of the matrix (i^) is CP violating. 
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with 

F(r) = ^(0, r, r) = ^o(r, r) = /' ^^V" ^^ . In ^^l^rfy . (20) 

Jo r-y{l-y) t 

Note that, as expected, if both M^^ = M^j for (g = t, 6) CP violation disappear in our 
approximation. This explains the cancelation in Eg . ([T9|) . 

Similar contribution from the Z exchange diagram can be derived based on the inter- 
action, 

Cz = ^^Z (h^iTf - Cf sin" dw)f + eUtilfi + ~}fJ\fbTlfh] . (21) 

cosOw \ ■< ' J 

j^, ^ ( Tl cos^ e, - e, sin^ 9^ \Tl sin 26, \ 

\ \Tl sin 29, Tl sin^ 9, - e, sin^ 9w ) ' 

3acm tan/5(|T/^ — ee sin^ 9w) me 



' ^02 327r^ sin^ 6^1^ cos^ 9w M\q 

Although we explicitly use the lepton i for our EDM calculation, the formalism can be 
easily extended to the quark. The color EDM of quark can also be obtained from Eq. (|T^) 
by changing color factors. After including the renormalization group evolution factor, we 
use the simple nonrelativistic quark model to estimate the neutron EDM0]. 

In Fig. 3 we show contributions to the EDM of the electron and the neutron. For 
EDMs of the electron and the quarks, contributions from the charged Higgs boson (plus 
W boson) and from the pseudoscalar neutral Higgs boson (plus either photon and Z 
boson) are all shown for comparison. For the EDM of the neutron, there is additional 
contribution due to chromoelectric dipole moments of quarks from two-loop diagrams 
with the pseudoscalar neutral Higgs boson plus gluons. For simplicity, we assume the 
/i parameter is real and Mj^ = M?^ = 0.6 TeV, At = Af, = i{l TeV) as in Ref.||. In 
Figs. 3a and 3b, we display the tan/5 dependence of each contribution to the EDM for the 
electron and the neutron respectively for both Ma = 150 GeV and 300 GeV (/i = 1 TeV). 
In Figs. 3c and 3d we display the fi dependence for the same choice of parameters with 
tan/5 = 20. In the minimal supersymmetric Standard Model (MSSM), the masses of the 
charged Higgs boson and the pseudoscalar boson are related by 

M^± =Mlo+M^ + e , 

where e = 3g2mf\ia{l + (m^/m^))/(87r^M^) is the the quantum correction^]. We have 
used the tree level result in plotting Figs. 3a-d. However, if quantum corrections are in- 
cluded, Mh± and Ma become less closely related. Therefore in Figs. 3e and 3f we display 



the dependence on Mh± and the dependence on M4 for their respective contributions. 
Overall, the EDMs are not very sensitive functions of Mh± and Ma- For the EDM of 
the neutron, the chromoelectric dipole moment contribution (due to pseudoscalar boson 
exchanges) still dominates over other sources. For the EDM of the electron, the contribu- 
tion of the charged Higgs boson exchange is in general small than that of the pseudoscalar 
exchange by about an order of magnitude. The A^Z exchange diagrams, which have not 
been included in Ref.0 contribute even smaller. 

Conclusion 

We have calculated the charged Higgs related two-loop contributions to the electric dipole 
moments of electron and neutron. For numerical simplicity, we ignore the generational 
mixing between squarks, however, it should be quite easy to incorporate if necessary. We 
find that the charged Higgs contribution are generally smaller than the neutral Higgs 
contribution calculated earlier^] in MSSM without the squark flavor mixing which we 
assume. However, one can imagine that in theories beyond MSSM, the two contributions 
may involve independent parameters that should be constrained separately. It is straight- 
forward to generalize our analytic result to accommodate theories beyond MSSM or to 
include squark flavor mixing. 



While we were preparing this manuscript, we became aware of a preprint |jTO| which 



aims to calculate the same contribution. Our analytic results in Eq. (|]) and Eq. (14) differ 



from those in Ref. ||10||. In particular our amplitude in Eq. (|T^) has the gauge invariant 
form, while the corresponding one in Ref. |TU[ is not. This difference does not affect the 
order of magnitude of the numerical results very much. However, the result for the charged 
Higgs contribution to EDM of electron in Ref. [|T0[] starts decreasing for large tan/5(> 35), 



but not in our gauge invaraiant result as shown in Fig. 3a. 
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Fig. 3a Numerical estimates of various contributions to the EDM of the electron as a func- 
tion of tan/3, for the case M?^ = M?^ = 0.6 TeV, Ma = 150 GeV and 300 GeV, 
real /U = 1 TeV and imaginary At = Ai, = i{l TeV). 



Fig. 3b Numerical estimates of various contributions to the EDM of the neutron as a func- 
tion of tan/3, for the case M?^ = M?^ = 0.6 TeV, Ma = 150 GeV and 300 GeV, 
real /^ = 1 TeV and imaginary At = Ab = i{l TeV). 
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Fig. 3c Numerical estimates of various contributions to the EDM of the electron as a func- 
tion of /i, for the case tan/3 = 20, M?^ = M|^ = 0.6 TeV, Ma = 150 GeV and 
300 GeV, real fi and imaginary At = Ah = i{l TeV). 

Fig. 3d Numerical estimates of various contributions to the EDM of the neutron as a func- 
tion of /i, for the case tan/3 = 20, M?^ = M|^ = 0.6 TeV, Ma = 150 GeV and 
300 GeV, real /i and imaginary At = Aj, = i{l TeV). 
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Fig. 3e Numerical estimates of H^ contributions to the EDM of the electron as a function 
of Mh±, for the case M? = M~ = 0.6 TeV, real /i = 1 TeV and imaginary 
At = At = 1(1 TeV). 



Fig. 3f Numerical estimates of A^ contributions to the EDM of the electron as a function 
of Ma, for the case M? = M? = 0.6 TeV, real /i = 1 TeV and imaginary At = 
Ab = i{l TeV). 
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